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The final state interaction (FSI) in Deep Inelastic Scattering (DIS) of leptons off 
a nucleus A, due to the propagation of the struck nucleon debris and its hadroniza- 
tion in the nuclear environment is considered. The effective cross section of such a 
partonic system with the nucleons of the medium and its time dependence are esti- 
mated, for different values of the Bjorken scaling variable, on the basis of a model 
which takes into account both the production of hadrons due to the breaking of the 
color string, which is formed after a quark is knocked out off a bound nucleon, as well 
as the production of hadrons originating from gluon radiation. It is shown that the 
interaction, the evolution and the hadronization of the partonic system in the nuclear 
environment can be thoroughly investigated by a new type of semi-inclusive process, 
denoted A(e,e'{A — 1))X, in which the scattered lepton is detected in coincidence 
with a heavy nuclear fragment, namely a nucleus A — 1 in low energy and momen- 
tum states. As a matter of fact, if FSI is disregarded, the momentum distribution 
of A — 1 is directly related to the momentum distribution of the nucleon before 7* 
absorption, i.e. the same quantity which appears in the conventional A{e, e'N)X 
process, where N denotes a nucleon. The rescattering of the struck nucleon debris 
with the medium damps and distorts the momentum distributions of A — 1 in a way 
which is very sensitive to the details of effective cross section of the debris with the 
nucleons of the medium. The total cross section of the process A(e, e'{A — 1))X on 
^He, ^^O, and '^^Ca, related to the probability that after a target nucleon experi- 



ences a DIS process, the recoiling A — 1 nucleus remains intact in spite of the strong 
FSI, is evaluated, and the distorted momentum distribution of the recoiling A — 1 
system is obtained. It is shown that both quantities are very sensitive to the details 
of the early stage of hadronization of the nucleon debris in the nuclear medium. 
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I. INTRODUCTION 



Lepton scattering off nuclei in the Deep Inelastic Scattering (DIS) regime represents a 
powerful tool to investigate a wide range of physics processes related to strong interaction 
physics which are more difficult to study in lepton scattering off free nucleons. As a matter 
of fact, in lepton nucleus scattering the interaction of the debris of the struck nucleon with 
the nuclear medium during the hadronization processes, could provide in principle valuable 
information on the space time structure of the hadron formation mechanism (see e.g. 
Among others, we should stress two main motivations for thoroughly investigating quark 
interaction in the hadronic medium: 

1. the first one, as already pointed out, is related to the possibility to understand the 
very mechanism of the formation time and hadronization; 

2. the second one stems from the obvious necessity, once a workable model for quark 
propagation and hadronization is developed, to apply to the treatment the Final State 
Interaction (FSI) in various processes involving nuclei, like e.g., ultra relativistic heavy 
ion collisions, aimed at observing a possible quark-gluon plasma formation, and semi 
inclusive lepton DIS scattering processes, aimed at investigating possible distortions 
of the nucleon structure function of a bound nucleon (EMC effect). 

It is the aim of this paper to update a comprehensive and workable model to treat the 
propagation and re-interaction in the medium of a nucleon debris produced in a DIS process 
off a bound nucleon [Q, and to apply it to a recently proposed semi-inclusive DIS process on 
nuclei |p, which could be useful not only to investigate in a more detailed way the mechanism 
of formation length and hadronization, but also to obtain more reliable information on EMC- 
type effects. 

To date, the information on hadron formation length comes mainly from the measurement 
of the multiplicity ratio of the lepto-produced hadrons in semi inclusive A{e, e'h)X processes 
([S H) D)) whose interpretation on the basis of the quark re-interaction model we consider in 
this paper appeared to be very convincing (see e.g. for a discussion of recent HERMES 
data). However, it should also be pointed out, that more exclusive processes, e.g. of the type 
we are going to discuss, though difficult to perform, could in principle provide more direct 
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information on quark re-interaction and hadronization mechanisms. As for the EMC-type 
effects, in spite of many experimental and theoretical efforts (for a recent review see P]), 
the origin of the nuclear EMC effect has not yet been fully clarified, and the problem as 
to whether, and to which extent, the quark distributions of nucleons undergo deformations 
due to the nuclear medium, remains open. This is why various semi-inclusive experiments 
in which another particle is detected in coincidence with the scattered electron have been 
proposed. Most of theoretical studies in this field concentrated on the process D{e, e'N)X, 
where D denotes the deuteron, a nucleon, and X the undetected hadronic state. Current 
theoretical models of this process are based upon the impulse approximation {I A) (or the 
the spectator model ), according to which: i) X results from DIS off one of the two nucleons 
in the deuteron, ii) the second nucleon recoils without interacting with X and is detected 
in coincidence with the scattered electron (for an exhaustive review see [^). The model has 
been improved by considering that the detected nucleon could also originate from quark 
hadronization |1T0| , |Tl[] , and has also been extended to complex nuclei by considering the 
process A{e, e'N)X, by assuming that DIS occurs on a nucleon of a correlated pair, with 



the second nucleon A^ recoiling and detected in coincidence with the scattered electron ^1 
In all above calculations, however, the nucleon debris created by the virtual photon is 

assumed to propagate without re-interacting with the spectator nucleus, which, therefore, 

always remains intact, an assumption which, at first sight, might appear unjustified. 

As a matter of fact, DIS scattering off a bound nucleon results in the production of a 

multi-particle final state with an effective mass squared equal to , 



s' -m%-Q' + 2mNiy-2^/iy^~Q^pL 

= Q\--l)+ml-2\ci\pL (1) 

X 

where = — is the four- momentum transfer, u the virtual photon energy in the rest 

frame of the nucleus, pi the longitudinal Fermi momentum of the nucleon relative to the 

Q2 

direction of the virtual photon (q 11 z), the nucleon mass, and x = the Bjorken 

scaling variable (we neglect here the binding energy of the nucleon). At high energies and 



far from the quasi-elastic region {x ~ 1), the effective mass is large, ^ m^v, and one 
could expect production of many particles which can interact traveling through the nucleus. 
This would substantially suppress the probability for the spectator nucleus to remain intact. 



However, the process of multi-particle production has a specific space-time development, and 
it turns out that not so many particles have a chance to be created inside the nucleus. 

Recently , a new type of semi-inclusive process on complex nuclei has been considered, 
namely the process A{e, e'{A — 1))X, in which DIS occurs on a mean-field, low-momentum 
nucleon, and the nucleus {A — 1), recoiling with low momentum and low excitation energy, is 
detected in coincidence with the scattered electron. Within the lA, it has been shown that 
such a process exhibits a series of very interesting features which could in principle provide 
useful insight on the nature and the relevance of quark FSI in DIS off nuclei, the validity of 
the spectator mechanism, and the medium induced modifications of the nucleon structure 
function. In the present paper we go beyond the lA by considering the effects of the quark 
re-interaction in order to clarify: i) to which extent the conclusions reached in Ref. |^ will 
be affected by the FSI, and ii) if, and to which extent, the process is sensitive to the details 
of quark hadronization in nuclear environment. 

Our paper is organized as follows: in Section 2 the formalism of FSI is presented; the 
application of the theory to the process A{e,e'{A — 1))X is illustrated in Section 3; the 
Summary and Conclusion is given in Section 4. 



II. FINAL STATE INTERACTION AND HADRONIZATION IN 
SEMI-INCLUSIVE PROCESSES 

A. Coherence time for particle production. The color string model 

After a quark is knocked out off a bound nucleon by the virtual photon, a color field is 
stretched between the quark and the remnants of the nucleon. In the color fiux tube model 



T2| , |T^, one assumes that this process is adiabatic, i.e. gluon radiation is neglected while 
the color field is squeezed by the QCD vacuum to a color tube of a constant cross area. 
Assuming that the transverse dimension of the tube is much less than its length one can 
call it color string. 

The important parameter of the model is the string tension, i.e. the energy density per 
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unit of length, which is related to the hadronic mass spectrum ||T2 



1 IGeV 



(2) 



27ra'^ fm 



where ~ 0.9GeV^~^ is the slope of the leading Regge trajectories. 

The energy of the leading quark degrades with a constant rate, dE/dz = —k, {z is 
the longitudinal coordinate) which is invariant relative Lorentz boosts. At the same time, 
the rest of the nucleon {e.g. a diquark) speeds up with the same acceleration, so the total 
momentum and energy of the debris of the nucleon remain constant. The string itself carries 
only energy, but no longitudinal momentum. 

Naively one might think about a long string stretched across the nucleus. However, the 
slow end of the string is accelerated and reaches soon the speed of the light. It turns out 
that the maximal length of the string in the rest frame of the target nucleon is. 



where m^g is the mass of the rest of the nucleon, which we conventionally call a diquark. 
Assuming this mass to be less than the mass of the nucleon, one gets L^ax < 1 Thus 
the string propagating through the nucleus is rather a short object in the nuclear rest frame. 
Moreover, subsequent decays of the string via spontaneous qq pair production from vacuum 
make the string even shorter (see below). 

Since the interaction cross section of a high-energy colorless object in QCD depends only 
on its transverse size, one can assume that for such a string it should be of the order of 
nucleon-nucleon total cross section. 

An important phenomenon related to the string evolution is the spontaneous creation of 
quark-antiquark pairs from vacuum via the Schwinger mechanism. Since the string potential 
is a linear function of the distance, a created qq pair completely screens the external fields 
in between, therefore, it breaks the string into two pieces. If so, the interaction cross section 
may nearly double compared to one string. This is important and must be taken into 
account. 

The probability W{t) for a string to create no quark pairs since its origin is given by 





6 



where w is the probabihty rate to create a hght qq pair within a unit of length of the string 
and a unit of time {t = z, in natural units), and L{t) is the time dependent length of the 
string. Note that Eq. (^) is invariant relative to longitudinal Lorentz boosts. 

The key parameter w can be estimated either using the Schwinger formula, or by cal- 
culating the decay width of heavy resonances |jl2|, [13|, both ways giving w ^ 2 fm~'^. One 
can also evaluate w from the momentum distribution of the recoil protons in the reaction 
pp —>■ pX |T5[. If the final proton stays in the fragmentation region of the target, it 



acquires momentum due to the acceleration by the string up to the moment At of the first 
string breaking. 

This moment is determined according to (^) by the condition 

(5) 

The momentum which the proton gets during this time interval, p = KAt, is related to the 
Feynman xp of the leading proton (in the anti-laboratory frame), p = m^v (1 — x\)/2xf- 
The mean value of xp is known from data, {xp) ~ 0.5. Therefore, {p) ~ IGeV/c, and 
w = ^ 2k^/p'^ ^ Thus, all of these estimates converge at about the same 

value. 

From this consideration we found that the mean time of breaking of the string after 
its production is At 1 fm. Since after each breaking the leading piece of the string 
becomes twice shorter (in the lab frame), the time interval up to the next breaking, and 
correspondingly the momentum of the next produced particle, double. This is the way 
how string decay produces sequences of hadrons whose momenta are ordered in geometrical 
progression, corresponding to a plateau in rapidity scale. This bunch of particles with 
multiplicity rising with time propagates through nuclear matter with increasing probability 
to interact, i.e. to break up the recoiling nucleus. The effective interaction cross section of 
the partonic system developing in nuclear matter rises as function of time as, 

^e//(t) = a,r + ail,f nM(t), (6) 

where 

, , ln(l + t/At) ,_, 

nuit) = , (7) 
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and we have assumed that the string at the early stage before the first breaking interacts 
with the nucleon cross section. Subsequent decays of the string are assumed to be, str — >■ 
B + str ^ B + M + str B -\- 2 M + str. . ., where B and M are baryon and meson 
respectively. All other produced particles are mesons which are assumed to interact with 
the pionic cross section, a^^/^ = cr^^ ■ We also assume that they decay predominantly outside 
the nucleus. The latter assumptions is well justified for few-GeV mesons. 

The effective cross section grows logarithmically with time and in a long time interval, 
t oc Eg, when the hadronization is completed, reaches the maximal value, 

ama. ^ aff + <f {um) , (8) 

where (tim) is the observed mean multiplicity of produced mesons. If the energy of the 
quark initiating the jet is sufficiently high, the late stage of string hadronization happens 
outside the nucleus. 



B. Gluon bremsstrahlung 

An intuitive QED analogy for the string model would be a capacitor whose plates are 
much larger than the distance between them. These plates moving apart are losing energy 
mainly for creation of the static electric field, while photon radiation from the edges is a 
small correction. This analogy gives a hint to the restrictions for the application of the 
string model. Probably, soft inelastic interactions can be treated this way, although the size 
of a constituent quark and the length of the color-flux tube ~ 1/m (up to the first break) 
are of the same order, and gluon radiation might be an important correction. However, a 
quark knocked out by a highly virtual photon in DIS has size ~ l/Q, much smaller than the 
length of the string. This is apparently a situation when the "edge effect" of gluon radiation 
plays major role 0. 

In order to estimate the effect of gluon radiation one can rely upon perturbative QCD 



methods 0, [16], 0. We are interested in the time-dependence of the amount of radiated 
gluons, therefore the coherence time of radiation is important. It depends on the quark 
energy which is about the energy of the incident virtual photon Eg ^ u , the transverse 
momentum of the gluon, and the fraction a of the quark light-cone momentum it carries 
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1,0, 

_ 2E,a{l-a) 

where is the quark mass which is not important for further estimates. 

The mean number of gluons which lose coherence and are radiated during the time 
interval t is given by TE\, 

Q2 1 

„„W=/d4 / ^o:^^B(.-U), (10) 

A2 kr/Eg 



where the number of radiated gluons as a function of a and kx reads 0, [19 

dnc 4as{k^) 1 



da dk^ 3 TT ak'^ 



(11) 



Here = An/d ln{k^/ Aqi^^) is the leading order QCD running coupling. We use the 

approximation of soft radiation, a <^ 1 and k"^ <^ . Since dnc/dk^ oc 1 / k"^ at k^ ^ Q"^ 
20| we can use as the ultra-violet cut off in the k^ integration in eq. ([To|). 



To avoid double counting we assume that at Q < A the string fragmentation mechanism 
dominates the nonperturbative dynamics of particle production. For this reason we intro- 
duced in (|10]) an infrared cut-off for the /c|,-integration, which should be taken of the 
order of the semi-hard scale characterizing gluon radiation. This scale was fixed in at 
A = 0.65 GeV from data on diffractive gluon radiation (the triple-Pomeron mechanism) in 
soft hadronic collisions. The relatively large value of A corresponding to a semi-hard scale, is 
dictated by the experimentally observed smallness of the cross section of single diffraction, 
pp — s> pX to states of large effective mass Mx- The parameter A controls the size of the 
virtual gluon clouds surrounding he valence quarks. Although the mean radius of the cloud 
ro — 0.3/m might seem to be small, it goes well along with the radius of gluon-gluon corre- 
lation, calculated on the lattice [^, or with the instanton phenomenology [^] which gives 
similar estimates. The appearance of such a semi-hard scale has presumably a nonpertur- 



bative origin and was interpreted in |T3 as a result of a strong nonperturbative light-cone 
potential between the parent quark and gluon. 



After integration in (10) we find that multiplicity of radiated gluons rises with time 
differently dependent on whether t is smaller or larger than to = 1/ {ttim xsj) = 0.2fm/xBj- 
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At t < to, 



ncit) 



16 



In 



Q 



+ In 



tA 



QCD 



27 { \X J V 2 
At it > to the t-dependence starts leveling off, 



In 



16 



+ In 




+ In 



tA, 



QCD 



In 



HQ/^qcd) 
ln(A/AQCD 

\n{Q/AQCD^/h^t) 
ln(A/AQCD 



(12) 



2 A, 



QCD 



In 



HQ/Aqcd 



(13) 



ln{Q/AQCDVh/t)_ 

and reaches the maximal constant value when the time exceeds the full hadronization time, 

t > ^oQVA^ 



n 



G 



16 
27 



-In 



In 



In 



MQ/Aqcd) 
\n{\/AQCD 



(14) 



,2 Aqcd 

Apparently, all the three regimes Eqs. (P!^)-([T^) match. 

Thus, we have arrived at a similar logarithmic growth of the amount of the produced 
particles with time as in the string model. Note that no gluons is radiated at Q ^ A, which 
is reasonable since at this value the onset of the nonperturbative dynamics of gluon radiation 
is expected Jl^. At Q > A the amount of gluons slowly rises with Q. This observation does 
not contradict the fact that gluons provide the main contribution to the energy loss at large 
values of • Although the amount of hard (high kr) gluons is small, each of them takes 
away a large value of energy , uj > k'^t. 

One can replace each radiated gluon by a color octet qq pair with accuracy and 
rearrange these pairs to form colorless qq dipoles. If to treat those dipoles as produced 
mesons, the effective cross section (0) get an extra contribution. 



a, 



tot 



+ a 



nN 
tot 



(15) 



which smoothly switches into the expression (|^) at Q X. The effective absorption cross 
section cie// as function of time is depicted in Fig. |l| for different values of xbj- 
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H — 1 — I — I — I — I — I — 1 — I — 1 — I 
2 4 6 8 10 

z, fm 

FIG. 1: The effective quark debris-nucleon cross section vs the longitudinal co- 
ordinate z in correspondence of = WGeV^ and two values of the Bjorken 
scaling variable x. The open dots corresponds to the string model (Eq.^), and 
the full dots and stars represents the total (Tg// i.e. the sum of the string and 
gluon radiation contributions (Eq. p!5|) . 

III. THE SEMI-INCLUSIVE PROCESS A(e,e'(A - 1))X 

A. The Impulse Approximation 

The process we are going to consider is the one in which 7*, interacting at high with a 
quark of a mean field nucleon (to be labeled by "1", and considered to be a proton), having 
four-momentum pi = (pio,Pi), with pi = — P^-i, creates a nucleon debris which propagates 
in the A — 1 nucleon system, which recoils with low momentum Pa-i and low excitation 
energy, and is detected in coincidence with the scattered electron. We are interested in the 
propagation and interaction of the quark debris with the nuclear environment, but, to better 
understand the problem, we first analyze the Impulse Approximation {I A), when any kind 
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of FSI is disregarded. In one-photon-exchange approximation, the differential cross section 
in the laboratory system has the following form 



dxdQ'^dPA-i 

= K\x, Q\ zf^)z'^^^F^I^{xAi Q\v\) P\E, \Pa-^). (16) 

where = —c^ = —{ke — k'^"^ = q"^ — = AEeE'^sin'^ ^ is the four-momentum transfer 
(with q = ke — k'^, V = Ee — E'^ and 9 = 0^,), x = Q'^/2Mv is the Bjorken scaling 
variable, pi = (pio,Pi)j with pi = —Pa-i, is the four momentum of the nucleon before 
interaction with 7*, F^^^ the DIS structure function of a nucleon bound in nucleus A, and 
K'^lx, Q^, uai 2^1^^) the following kinematic factor 



X \yA 



;i7) 



with 



y = v/E^ , yA = (pi ■ q)l (pi ■ k^) (18) 

- (^), z^ - Mu ' ^ ' 

^1 

In Eq. |TB|, the quantity P{E, I-P4-1I) is the Nucleon Spectral Function 

P\E, \Pa-,\) = Y.\{Pa-i, ¥a-M)\'KE - {Er^^n + E{_,)) (20) 

/ 

where and '^^a-i wave functions of the target nucleus and the final nucleus {A — 1) 

in excited intrinsic state /, respectively, E = Emm + E-^_-^ is the nucleon removal energy, 
i.e. the energy required to remove a nucleon from the target, leaving the A — 1 nucleus with 
excitation energy E^_^, and, eventually, Emin = M + Ma-i — Ma- 

From now-on, the overlap appearing in Eq. (pO]) will be called the transition form factor 
of the process and will be denoted as follows 

F{_,^^{Pa-i) ^ {Pa-i.K-^'^'a) 

„ A A 

= / e^''^-^^^^Li(^"'2---^"'A)<(n,rl...f^)5(^f,)]^rfr-;, (21) 

j=i i=i 
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Because of energy conservation 



u + Ma = P(a-i)o + Pxo (22) 

where pxo is the total energy of the debris, the removal energy can be written in the following 

way 

E = u + M-p,o (23) 

if the total energy of the system A — 1 is approximated by its non-relativistic expression 
and the recoil energy disregarded. 

As is well known |2^, the integral of the spectral function over the removal energy E 
defines the (undistorted) momentum distributions 



.^(i^)^ye--)p(,n^^.' 

= J dEP\E, IpI) = + <(p-) 



(24) 



where 



A 

p{r, r') = / ^Yir, r 2 . . . fA)^'^(f ', fa . . . f^) JJ ' (25) 

i=2 

is the one-body mixed density matrix. In Eq. ^ Uq represents the mean-field uncorrelated 
momentum distribution arising from the summation over the discrete hole states of the final 
system, whereas is the correlated momentum distribution resulting from the summa- 
tion over high excitation states of the final system, which originate from nucleon-nucleon 
correlations. 

In this paper we will consider semi-inclusive processes, when the cross section (^) is 
integrated over the removal energy E, at fixed value of Pa-i- Thus, owing to 

A 

^{-i\r, ' . . . r4 'WA_,{f2 . . . r-k) = n ^(^^ " ') ' (26) 

the cross section ( [T6| ) becomes directly proportional to the momentum distribution 
n(|P4„i|)^. However, since, as previously stated, we will only consider the formation of 
a low momentum and low excitation energy A — 1 system, the summation over / is effec- 
tively limited to the hole states a of the initial nucleus, which means that the only relevant 
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quantity is the low-momentum part of 

^o(i^A-ii) = $^m</i(iPA-ii)r , (27) 

a<F 

where the sum extends over the states below the Fermi sea occupied in the ground state. 
Following 0, let us briefly discuss the A dependence of the cross section (|T^). In this 
respect it should be pointed out that nuclear effects are not only generated by the nucleon 
momentum distribution ?7,^(|P4_i|), but also by the quantities ha and z['^\ which differ from 
the corresponding quantities for a free nucleon {y = v/ E^, and z^'^ = 1) if the off mass shell 
of the nucleon {p\ ^ ) generated by nuclear binding is taken into account. However 
such a dependence is not only very small, but it completely disappears in the ratio between 
the cross sections from nuclei A and A' , 

RbA^Bj/Z^ ,Q ,\Fa-iIA,A) - / {A') A>(\?, ^ 

F2 i.XBj/z\ \Q^)n^ {\Pa-i\) 



<i\PA- 



I), (29) 



<i\PA-l\) 

which means that in the Bjorken limit the A dependence of the ratio R is entirely governed 
by the A dependence of the nucleon momentum distribution r;,^(|P4„i|) in nuclei A and A' . 
Since, as shown in Fig. ^ exhibits a strong A dependence for low values of I-P4-1I, a 
plot of R versus |-P4_i| should follow the behavior of nQ{\PA-i\) in nuclei A and A' and the 
experimental observation of such a behavior would represent a stringent test of the spectator 



N/A 
2 



mechanism independently of the model for 

The expected behavior of the ratio Eq. ( p9D for A = 2 and different values of A' is 
presented in Fig. ^ . These results clearly show that the observation of recoil nuclei in 
the ground state, with a |P4-i|-dependence similar to the one predicted by the momentum 
distributions, would represent, on one hand, an indication that the FSI between the lepto- 
produced hadronic states and the nuclear medium is such as to leave intact part of the the 
final {A — 1) nuclei. One could hope that the number of the detected nuclei together with 
the momentum dependence of the cross section could provide information on the nature and 
the details of the hadronization mechanism in a more sensitive way than e. g. the energy 
transfer dependence of forward hadro-production in the process A{e,e'h)X 
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10 -1 ■ 1 ■ 1 ■ 1 ■ 1 ■ r- 

0.0 0.2 0.4 0.6 0.8 1.0 

P,.,,fm' 

FIG. 2: The proton momentum distributions in (full), ^He (dotted), ^^O 
(dashed) and ^^Ca (dot-dashed) calculated using realistic nucleon-nucleon in- 
teractions (see for original references). 

B. Final state interaction and hadronization 

In this Section we are going to consider the modifications of the cross section induced by 
the FSI of the nucleon debris produced in the DIS on a bound nucleon. In this case, the 
Spectral Function has to be replaced by the Distorted (D) Spectral Function, which can be 
written in the following way 

P^iE,P^_,) = J2\FA-UPA-i)mE-iE^,^ + E{_,)) , (30) 
/ 

where the transition form factor is 




Slin...rA)'^{^i\r2-.-rA)'^Airi,r2--.rA) 5 H^fi (31) 

Vj=l / i=l 
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10^ 



■ \ \ 



10' -I ■ 1 ' 1 ' 1 ■ r 

0.0 0.2 0.4 0.6 0.8 1.0 

FIG. 3: The ratio Ra,a'{Pa-i) (Eq. (|2|) with Pa-i = \Pa~i\) corresponding 
to ^ = 2 and ^' = 4 (dotted), A' = 16 (dashed) and A' = 40 (dot-dashed). 



The quantity So in Eq. (|3lD is the Glauber operator, which describes the FSI of the debris 
from the struck proton with the {A — 1) system, i. e. 

A 

Scin . . . f^) = n [l - r^*'^(^i - b^, - z^)Q{zi - zi)] , (32) 



i=2 



where hi and Zi are the transverse and longitudinal components, respectively, of the coor- 
dinate of nucleon "i", rj = {bi,Zi), r'^*'^{b) is the Glauber profile function describing the 
elastic scattering of the debris, denoted A^*, with the nucleons of the {A — 1) system, and 
the function 6{zi — zi) takes care of the fact that debris of the struck proton "1" propagates 
along a straight-path trajectory, so that they interacts with nucleon "j" only if Zj > Zi. 
We have chosen the longitudinal axis z along the momentum of the virtual photon, and 
besides the usual dependence on the transverse separation bi — bi between A^* and Ni for a 
high-energy amplitude, we have also reserved a dependence of T^*^ on the longitudinal sep- 
aration Zi — Zi, which should take care of the time dependence of the effective cross section 
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discussed above. Let us eventually point out that because of the effects from the FSI, the 
Distorted Spectral Function depends now upon the vector Pa-i 7^ Pi- With the FSI taken 
into account, the Distorted cross section is now Eq. [l^ with P^{\Pa-i\, E) replaced by 
P^''^(P4_i, E). The cross section integrated over the removal energy becomes proportional 
to the distorted momentum distribution 



n 



D 



A-i) = E|^i;A-i(^^-i) 



(2vr)- 



-3 / APA-i{r-f") 



Pd(^, r')drdr' 



(33) 



where 



A 

/ ^2 • • • rA)Sl{f, r2 . . . fA)SG{f', fa . . . fA)K{r", fa . . . f^) JJ dvi (34) 

i=2 



is the one-body distorted mixed density matrix. 

A general approach to calculate PD(f;f') has been developed in Ref. ||2^ in terms of 
correlated wave functions. Since we are interested in the low momentum part of the wave 
function, and also due to the exploratory nature of our calculations, we limit ourself here 
to a traditional Glauber-type nuclear structure approach in which \E'^(f, fa . . . ta) is written 
as a product of a function 0, describing the motion of nucleon "1" and the wave function 

A-iir2---rA) 



■ta) of the spectator, and \'^^j^_i{r2 ■ ■ ■ ^a)\'^ is factorized into a product of single 



particle densities. One obtains, in 


this case. 




n^{pA-i) 


= N{Pa-i) = 


FA,A-liPA-l) 


with 






F?,a-i(Pa-i) ^ ^ 


1 e^-^-^-i"0(f) 


' S{b,z) ' 
2{A-1)_ 


and 


00 




S{b, z) 


= j dz' pA(b,z 

z 


')aeff{z' -z) 



(35) 



A-l 



df 



(36) 



(37) 



where pA{b,z) is the nuclear density (normalized as J d^r pA{f) = A), (j){f) describes the 
relative motion of the struck nucleon with respect to the spectator A — 1 nucleus, and 
(Jeff{z' - z) is given by (|) or (||). 
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We have estimated the Distorted form factor, Eq. (^7]), within the optical approximation, 
obtaining 



oo 

1 



(38) 



where Pt and Pl denote the transverse and longitudinal components of Pa-i with respect 
to q. The optical approximation is a rather good one for heavy nuclei, but it should be 
abandoned in case of light nuclei for which Eq. (0) has to be used. 

The total cross section of the process is proportional to the integral of the distorted 
transition matrix element, obtaining 

oo 

(Ttoioc J (fPA-i \fI^_^{Pa-i)\^ = j d% J dzp{b,z)exp[-S{b,z)'\ . (39) 

— oo 

IV. RESULTS OF THE CALCULATIONS OF THE TOTAL CROSS SECTION 
AND THE DISTORTED MOMENTUM DISTRIBUTIONS 

We have calculated both the total cross section and the distorted momentum distribution 
of the process A{e, e'{A~l))X off ^ife, ^^O, and ^^Ca; the hit nucleon is considered to be a 
proton. For the nuclear density pA{b, z) we have used both the Harmonic Oscillator and the 
Fermi distribution forms. As for the function (f){bi,zi), which describes the relative motion 
of nucleon 1 and the spectator nucleus A — 1, it has been chosen such that when FSI are 
absent, A^(|Pa-i|) coincides with the low momentum part of the distributions given in Ref. 



| 25| . The parameters used for the calculation of the effective cross sections [Eqs. (|]) and 
dH)] were as follows: aff = 40 mb, a^^ = 20 mb, Aqcd = 0.25 GeV, X = 0.65 GeV. The 
results of the calculations are presented in Figs. |^-|^. 

Fig. ^ shows the ratio of the distorted to undistorted total cross sections defined by 
Eq. ( PSP calculated either taking into account {S{b,z) 7^ 0), or disregarding {S{b,z) = 0) 
the effect of the nucleon debris rescattering in the medium, respectively. 



As anticipated, the rescattering effects exhibit a decreasing A-dependence. The effect 
of gluon radiation, responsible for the difference between closed and open dots, is a nearly 
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FIG. 4: The ratio between the total cross section (Eq. ( |39| ) calculated tak- 
ing into account {S{b,z) ^ 0, otot = (^tot^) disregarding (5(6,2;) = 0, 
atot = '^tot^) Final State Interaction. The open dots correspond to the 
effective cross section of the nucleon debris given by the color string model 
[Eq. (^)], whereas the full dots correspond to the cross section where the 
gluon bremsstrahlung has also been considered [Eq. (|l|)] at = lOOGeF^. 
The stars represent the proton transparency calculated in Refs. p^ , [2^ for the 
reaction A{e^ e'p)X. 



A-independent 10% correction. We have also calculated the nuclear transparency for the 
reaction of quasi-elastic scattering A{e, e'p)B where the intact struck proton (Glauber 
approximation) propagates through the nucleus. These results depicted by star points 
clearly exhibits the expected weaker absorption effects for the proton compared to the 
debris in the process A{e, e'{A — 1))X. The momentum distributions, Eq. (pTf), for "^He, ^^O 
and ^°Ca, are compared in Figs. ^|-|^ with the distorted momentum distributions defined 
by Eq. ( |35D calculated for Pt = 0. 
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FIG. 5: The proton Momentum Distribution for '^He shown in Fig.l (NO 
FSI) compared with the Distorted Momentum Distribution N{Pa~i) (Eq.(|35|)) 
plotted vs Pl for Pt = 0. The curve labeled by open dots has been obtained 
using the effective cross section for the nucleon debris corresponding to the 
color string model (Eq. @), whereas the other curves correspond to the cross 



section which includes also the gluon bremsstrahlung (Eq.(15)). The stars 
represent the distorted proton momentum distributions calculated in Ref. [^] 
for the semi-inclusive process ^He{e, e'p)X . 



Concerning the results presented in these Figures, the following remarks are in order: 

1. Due to the smaller dimensions, rescattering effects are less important in '^He than 
in heavier nuclei. As a matter of fact, in the former case, they simply reduce, in 
the considered range of momenta, to an almost constant reduction of about 50 %, 
whereas it appears that for medium weight and heavy nuclei not only the absorption 
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FIG. 6: The same as in Fig.^ but for ^^O. The proton rescattering in the 
reaction ^^0{e,e'p)X damps the momentum distribution (fuh curve) by a 
factor of about 0.5 



is substantially stronger, but appreciable distortion effects are apparent, with the role 
of gluon radiation increasing with A. 

2. It is interesting to compare the process A{e, e'{A — l))X we are investigating occurring 
at Bjorken x <^ 1, with the semi-inclusive nucleon knock-out A{e, e'p)X occurring at 
Bjorken x ~ 1 (note that X refers to the proton debris in the A{e, e'{A — l))X process, 
and to a y4 — 1 nucleon state in the A{e,e'p)X process). Within the PW approxi- 
mation the two processes are proportional to the same nuclear part, viz the nucleon 
momentum distribution n'^(p) (Eq. (0)), whereas when the FSI is considered the first 
process will be distorted by the nucleon nucleon debris rescattering, and the second one 
by the proton rescattering. Thus, it appears that by comparing the two processes the 
differences between nucleon debris and proton propagations in the nuclear medium can 
be investigated. The calculation of rescattering effects in the processes ^ife(e, e'p)X, 
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FIG. 7: The same as in Fig. ^ but for "^^Ca. The proton rescattering in the 
reaction ^''Ca(e, e'p)X damps the momentum distribution (full curve) by a 
factor of about 0.5 [|^. . 
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0(e, e'p)X and ^'^Ca(e, e'p)X has been performed in by a Glauber-type approach 
and the results can be summarized as follows: in the momentum region considered in 
this paper, the effects of the FSl, due to the proton rescattering, damp the momentum 
distributions by an almost momentum-independent amount, the reduction factors be- 
ing 0.85, 0.5 and 0.5 in ^ife, ^^O and ^^Ca, respectively. By comparing these results 
with the ones presented in Figs. we can conclude that, apart from ^He, the effects 
of proton and nucleon debris propagation expected, very different. 

In Fig. H the ratio Ra,a' (Eq- (HD)) which is plotted in Fig. ^ in the case of the PW 
approximation, is shown when the nucleon debris rescattering is taken into account in the 
final state. It is gratifying to see that when the final state rescattering of the debris is 
taken into account, the differences between different nuclei, i.e. nuclear effects, are even 
emphasized. 
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FIG. 8: The same quantity as in Fig. |3| obtained with the Distorted Momen- 
tum Distributions shown in Figs.^-^ (STRING). Dotted hne: ^He; dashed 
hne: ^^0\ dot-dashed hne: '^^Ca. 



In Fig. ^ the dependence of our results upon the choice of the function (j){b, z) describing 
in the target ground state the relative motion between the hit nucleon and the spectator 
y4 — 1 , is exhibited by comparing the results obtained with the realistic function used in 
the calculations with the results obtained with a simple Gaussian function. It can be seen 
that, apart from ^i/e, the use of a Gaussian function is not recommended. 

Eventually, in Fig. |10|, the sensitivity of our results upon the quantities entering the 
effective cross section (Eq. |15D, viz a^l' , and A, is illustrated. 



V. SUMMARY AND CONCLUSIONS 



In this paper we have considered a theoretical model to describe the FSI arising from 
the propagation and hadronization of the partonic debris produced from the hard scattering 
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FIG. 9: The distorted momentum distributions for '^He (full), ^^O (long 
dashes) and ^'^Ca (short dashes), calculated with two different forms for the 
function (j){b, z) which describe in the target nucleus the relative motion be- 
tween the hit nucleon and the spectator nucleus (A— 1) (cf. Eq. (^). Gaussian 
form: open dots; realistic form used in this paper: stars. 

of a lepton off a nucleoli bound in a nucleus A. The FSI arises due to the rescattering 
of the hadrons which are formed both from string breaking and from gluon radiations. 
In order to experimentally investigate the correctness of the model, we have analyzed the 
semi-inclusive process in which, instead of a secondary, e.g. leading hadron, the whole A — 1 
nucleus in low momentum and energy states, is detected in coincidence with the scattered 
lepton. In absence of any FSI, the momentum distributions of the A — 1 nucleus would be 
nothing but the momentum distributions of the hit nucleon and its integral will give the total 
number of surviving A — 1 nuclei; when the FSI is switched on, the surviving probability 
will be reduced and the momentum distributions will be distorted. We have shown that 
both effects sensibly depend upon the details of the rescattering and hadronization of the 
nucleon debris. By comparing the results for the process A{e,e'{A — 1))X with the much 
investigated process of nucleon knock-out, A{e, e'p)X, the differences between the nucleon 
and its debris propagations in the nuclear environment can be investigated. For example, we 
have found that whereas 60% of the knocked out protons off ^^O at = GeV^ escape the 
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FIG. 10: The sensitivity of the distorted momentum distributions of ^^O upon 
the values of the Nucleon-Nucleon and Pion-Nucleon elementary cross sections. 
The full curve represents the Plane Wave approximation, whereas the other 



curves were obtained using the debris effective cross section given by Eq.l5 . 



nucleus without interaction, for a debris of a struck nucleon the survival probability reduces 
down to 25%. Furthermore, the momentum distributions of the recoiling A — 1 nuclei turn 
out to be very sensitive to the details of the model for hadronization in nuclear environment. 
Up to now, information on hadron formation time and hadronization, has been collected 
mainly from the measurement of the ratio between the semi-inclusive cross section for leading 
hadron production for a nucleus A to that for the deuteron. This ratio increases with energy 
approaching 1 when the formation time of leading hadrons exceeds the nuclear size. In the 
process we are proposing, the effects of FSI do not vanish with energy, and vary from a 
factor of 2-3 in a light nucleus like '^He, to orders of magnitudes in heavy nuclei. Moreover, 
they generate peculiar and strong distortions of the nucleon momentum distributions. Thus, 
from a theoretical point of view the semi-inclusive process A{e,e'{A — 1))X appears to be 



25 



very promising, for nuclear effects not only manifest themselves in a relevant quantitative 
way, but also produce peculiar qualitative effects. 

An essential advantage of the process under discussion, compared to inclusive hadron 
production, is the possibility to study the early stage of hadronization at short formation 
times without being affected by cascading processes. Indeed, no cascading is possible if the 
recoil nucleus {A — 1) survives. At the same time, most of hadrons with small momentum 
produced in inclusive process A{e, e'h)X originate from cascading of more energetic particles. 
In order to analyze data and extract information on the early stage of hadronization, one 
needs a realistic model for cascading what is barely possible. 
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